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Abstract 

A method of deriving bounds on the weak meson form factors, based on 
perturbative QCD, analyticity and unitarity, is generahzed in order to fuhy 
exploit heavy quark spin symmetry in the ground state (L = 0) doublet of 
pseudoscalar (B) and vector {B*) mesons. All the relevant form factors of 
these mesons are taken into account in the unitarity sum. They are treated as 
independent functions along the timelike axis, being related by spin symmetry 
only near the zero recoil point. Heavy quark vacuum polarisation up to three 
loops in perturbative QCD and the experimental cross sections a{e~^e~ T) 
are used as input. We obtain bounds on the charge radius of the elastic form 
factor of the B meson, which considerably improve previous results derived 
in the same framework. 



1 Introduction 



Bounds on the charge radius of the elastic form factor of the i?-meson were recently 
derived in a number of papers The interest in this form factor comes from 

the fact that it coincides, in the large quark mass limit, with the renormalized Isgur- 
Wise function of the heavy quark effective theory The short distance and 

finite mass corrections are in this case much smaller than for the flavor changing 
currents involved in the semileptonic decays of the B meson into D or D*. Therefore, 
rigorous bounds on this form factor are of interest for testing various nonperturbative 
techniques applied for the calculation of the Isgur-Wise function. 

The method applied in Refs. based on previous works ^ starts by ex- 

ploiting the same input as the standard QCD sum rules, i.e. the QCD euclidian 
expansion of a polarization function, related by analyticity and unitarity to the 
physical states of interest. However, while in the usual formulation of the QCD 
sum rules one tries, by suitable methods, to enhance the contribution of the low 
energies in the dispersion integral and saturates the unitarity sum by the lowest 
lying resonances, in the approach proposed in the dispersion relation is written 
as a rigorous integral inequality for the modulus squared of the form factors of the 
physical states along the time like region. By using in addition the analyticity prop- 
erties of the form factors, this inequality is shown to constrain the behaviour of these 
functions or their derivatives near the zero recoil or other points of physical interest. 

In refs. |T|-[^ the method was applied to the elastic form factor of the pseu- 
doscalar B meson. An attempt to exploit heavy quark symmetry in the ground B 
meson state doublet was made in [|l|, where the BB and BB* + B*B intermediate 
states were included in the unitarity sum, with the additional assumption that the 
relevant form factors of the B and B* mesons are identical along the whole unitarity 
cut. However, this is an unjustified extension of the heavy quark spin symmetry, 
which is valid only near the zero recoil point. As illustrated in [§| by specific models, 
the B and B* form factors can be indeed quite different along the time like axis, 
especially near thresholds. The problem was correctly solved in [^], where by means 
of special techniques allowing the simultaneous treatment of several analytic func- 
tions [T^,[jlT|, the inclusion the form factor of the B ^ B* transition was possible 
within the strict heavy quark spin symmetry hypotheses. More precisely, the elastic 
form factor and the BB* form factor were treated as distinct functions along the 
unitarity cut, being assumed to coincide only near the zero recoil point. This led to 
a considerable improvement of the bounds on the charge radius of the B elastic form 
factor: the range —4.5 < < 6.1, obtained in [Q] without imposing spin symmetry, 
was narrowed in |^ to —0.90 < < 2.60. 

However, in spin symmetry was not fully exploited, as the contribution of 
the B*B* intermediate states in the unitarity sum was not included. This problem 
is addressed in the present paper, where we treat simultaneously all the weak form 
factors of the B and B* mesons. The quadratic expression yielded by unitarity is 
written in a suitable "diagonal" form, which allows us to apply the optimization 
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theory for vector- valued analytic functions |T^ |TT| . The different thresholds in the 
unitarity sum and the subthresholds singularities of the various form factors are 
taken into account correctly. In this way the consequences of spin symmetry in 
the ground state doublet of the B mesons are exploited in an optimal way. The 
present paper contains in addition two important improvements of the work done 
before: we use as input the heavy quark vacuum polarization function computed 
in perturbative QCD up to three loops and we include in the dispersion 



relation for the polarization function the three T resonances with masses below the 
threshold for BB production (these terms were neglected in previous works 

In the next section we present the derivation of the bounds. Section 3 contains 
the numerical results and our conclusions. 

2 The derivation of the bounds 

We use the standard definitions of the form factors of the B and B* mesons 

< B{p')\V''\B{p) >= {p + pYF{q^) (1) 

< B*{p',e)\V^\B{p) >= e,p'^pf,V{q') (2) 



< B*{p', e')\V^\B*{p, e) >= F,{q'){e ■ e')P, + F,{q')[e,{e' ■ P) + e'^{e ■ P)] 
+Fs{q') ^''^^^f'^^ P, + F,{q')[eUe' ■ P) - e'^{e ■ p)]M!^Mfi , 



where = b'-f^b, e(e') denote the polarization vectors of the B* mesons, P = p + p' 
and q = p — p' . 

The form factors defined above have cuts in the complex plane t = q"^, from the 
threshold to ioi BB production to infinity. The effect of the lower branch cuts due 
to light intermediate states {7i'7i,KK, etc) is negligible [||. The three resonances 
T(IS'), T(2S'), T{3S) with masses lower than 2mB produce additional singularities, 
which can be approximated by poles on the real axis below to 0- On the other 
hand, heavy quark symmetry predicts definite relations among the form factors (|l|)- 

near the zero recoil point w = 1 {w = v ■ v', v and v' being the velocities of 
the initial and final meson, respectively). In this region some of the form factors in 

and (|^) are approximately equal to the elastic form factor (|T]) and other vanish. 
Specifically, for w ^ 1 one has 

V{w) = -FiH = F2H = F{w) , F3H = F4H = , (4) 

and we recall that F{w) satisfies the normalization condition 

F(l) = l. (5) 
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We are interested in finding restrictions on the slope of tliis function at zero recoil, 
or the so called charge radius, defined as 



= -F'{1\ (6) 

which differs by ||logas(m{,) [1^ from the charge radius p^^r of the universal Isgur- 
Wise function |]^. 

As in the derivation of the usual QCD sum sules, for studying the form factors 
we start by considering the vacuum polarization tensor due to the current V^: 

W^[q) = [q^^q- - g^'y^)I{{q^) = i J dxe'"" < 0\T{V^'{x)V'' {0))\0 > . (7) 
The first derivative of the invariant amplitude n(g^) satisfies the dispersion relation 

uv)-\r^dt, (8) 

the spectral function being defined by the unitarity relation 

{q'^q^ - g^''q')lmUit + te) =\T. j dpri2n)'5^'\q - pr) 

X < o|v''(o)|r >< r|v'^(o)+|o > . (9) 

Here the summation is over all possible hadron states F with appropriate flavor 
quantum numbers, with an integral over the phase space allowed to each interme- 
diate state. We shall include in this sum the three T resonances with masses lower 
than the threshold of the BB production and the contribution of the two-particle 
states \BB >, \BB* + B*B > and \B*B* > above this threshold (the T(4S') res- 
onance is not included, in order to avoid double counting 0). This contribution 
can be evaluated in a straightforward way, by using the definitions (0-0) of the form 
factors, performing the phase-space integration and the summation over the polar- 
izations of the B* intermediate states. Taking into account the positivity of the 
spectral function of 11, which follows from (^) we obtain the following inequality: 

1 27 

-Imn(t + ie) > y Mr-Tr-S(t - Ml ) 

n 



^ { (i - f)'' \mfQ{t - to) + (1 - - 'jf^'fjvmit - Q 

e{t-t*A , (10) 



+ (1 - ^)3/2 



At ^ Af 

mit)? + TrAm? + + (T;^)'|i^4(t)P 

In In 



where 

m = i^- mit) + ^F,it) + - mit). (n) 
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In (^) the widths Ft^ are defined through the parametrization 



T 



(12) 



of the cross section for T production, to = 4m^, = (m^ + mB*Y tg* = 4mg* 
are the thresholds for BB, BB* and B*B* production, respectively. We used the 
notation ti = {rriB* — ^-nd = 3 is the number of light quark fiavors which 

give identical contribution in the unitarity sum 0. It was convenient to write the 
contribution of the last four form factors in a "diagonal" form, as a sum of moduli 
squared of functions with the same analyticity properties as the original form factors, 
which will allow the application of the mathematical technique presented below. 
This calculation was rather tedious and for simplicity these terms were omitted in 
the previous paper devoted to this problem. 

By combining the dispersion relation (^ with the unitarity inequality (p!0| ) we 
obtain the following integral condition for the form factors of interest: 



n'(g') > 



487r2 



it-q 



2\2 



1 



I) 



3/2 



\F{t)\^dt 



+ 



(t-g2)2^ t 



1 - -^fiHi 



f* 9+ 

nx3/2f^|T/m|2, 



t 



-\V{t)\'dt + 



it-q 



2^2 



^0 ',3/2 
t ' 



Af - Af 

2|Fi(t)p + -\F,{t)\' + \F,{t)\' + {-fWm' 
^0 f-o 



where 



n'(g^) = U\q') - 
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ATTa" 



E 



iq' 



Mir 



dt\ , (13) 



(14) 



In the euclidian region < the function 11' (g^) can be calculated by applying 
QCD perturbation theory, with nonperturbative corrections included by means of 
OPE. Due to the large value of m^, the QCD expression of n'(g^) can be used also 
at = or even at positive much less than 4m^. Moreover, in this case the 



nonperturbative corrections are shown to be entirely negligible [Q. In the previous 
works |Il|-0 only the lowest order (one-loop) perturbative polarization function was 
used as input in eq.(|l^) (the terms containing the T poles being also omitted). 
In the present analysis we introduce explicitely in (|14D the contribution of the T 
resonances, using the experimental information on Ty- . In the same time we use as 
input the expression of the polarization function up to three loops |ll2[-||14||: 



n'(g2) = n'(o)(g2 



+ ^£(/f!)n'(i)( 



+ 



TC 



(15) 



with the MS coupling as(/i^) defined in the conventional way. We use the standard 
expressions 



n'(o)(g2 



1 v{3-v^) 
Jo [l — q^x/Aml 



2\2 



(16) 



4 



n'(i)( 



f (3 — V 



2\ 



247rm^ Jo (1 — q^x/Aml 



' TT 


i; + 3 


i- - -)\ 


2v ~ 


4 


\2 47iJ_ 



dx, (17) 



with V = a/1 — X. As concerns the last term in (15), we shall use the Taylor series 
around = 0: 

the coefficients Cn being given in eq.(ll) of ref. Jl^ (we recall that they depend 
explicitely on the normalization scale fi). 

With the above expressions, the input entering (|13|) is completely specified and 
this inequality can be viewed as an integral quadratic condition for the form factors 
of interest along the unitarity cut. By applying standard techniques of analytic 



functions |T0|, extended to " vector- valued functions" (see |ri[] and the references 
therein) we shall obtain from this condition a quadratic inequality relating the values 
of the form factors and their derivatives at the zero recoil point. Using then the 
relations (§-|^) we shall finally express the derived inequality as a constraint on the 
charge radius (^. 

We first conformally map the cut t = plane onto the unit disk in the complex 
plane z, such that the unitarity cut becomes the boundary \z\ = 1. Actually, 
since the integrals appearing in ([13D have different thresholds, we shall use for them 
different conformal mappings . More precisely, we take 

z{t) = ^^"^ (19) 
v^o — t + v^o 



for the first integral in the right hand side of (13) and similar expressions, with to 
replaced by and tg*, respectively, for the second and the third integral. By the 
mapping (p!9D the threshold to becomes z = —1 and the zero recoil point w = 1 
(equivalent to t = 0, since w = 1 — is applied onto the origin z = 0. Similarly, 
using the mappings suitable for the other integrals in (^) as explained above, the 
thresholds and t^* become also z = —1 and the corresponding zero recoil point 
is applied on the origin. It is easy to see that the conformal mappings used for the 
second and the third integrals transform the threshold to into a point situated inside 
the unit circle, close to —1. By performing the above changes of variable, all the 



integrals in (p!3D become integrals along the same contour, i.e. the boundary z = e 
of the unit disk. 

It is convenient to introduce a compact notation by defining the following func- 
tions of the variable z: 

Mz)=F{t), f,{z)=V{t), 
Uz) = Fi(t), U{z) = F,{t), h{z) = F,{t), h{z) = F,{t), (20) 

where F^ is defined in (|ll]). Using the conformal mappings defined above, the 
normalization condition (§) and the definition (|]) of the charge radius, one can 
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show easily that the functions fi{z) satisfy the relations 

/,(0) = 1, /;(0) = -8p\ t = 1, ..5; /6(0) = 0, (21) 

the derivative being with respect with z. Moreover, following the standard technique 
presented in we shall define a set of functions analytic and without 

zeros in the unit disk, whose moduli squared on the boundary are proportional 
to the positive weights appearing in the integrals (|13D , multiplied by the Jacobian 
1^1 of the conformal mapping These functions can be constructed easily and 
unambigously, by using the relations 

{i-zr \ t) ~ 8 ' 

dt 1 + ^ 1 r ^ (22-) 

dz ° (1-^)3' (t-g2) \l-d) to{l-zdy' ^ ' 

which follow from (p!9D, with 



v^o - + v^o 

With these defintions, we can write the inequality ([TsD in the equivalent form 

^ rj:\uo)hm'de<i (24) 

27: Jo ^ 



where the functions (jyiiz), obtained using ( ^21) can be written in a compact form as 
The parameters entering this expression are as follows: 



MO) = «i = 2, &i = 1/2, ci = 2, = ci 

02(0) = 2V20i(O), as = 2, 62 = -3/2, ca = 2, ^2 = 

= i^/sS- = = = * = 

04(0) = 2\/203(O), 04 = 2, 64 = -3/2, C4 = 2, 4 = rf* 

05(0) = as = 2, 65 = 1/2, C5 = 2, 4 = c?* 

06(0) = 8x/203(O), as = 2, foe = -3/2, cg = 2, d^ = d* , (26) 

with n' defined in (|14D and c? in ( ^3]) (d* is obtained from d by replacing by 
rriB*). 
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As discussed above, the form factors appearing in (|T3|) have three poles on the real 
axis below the threshold to = 4:m%, due to the three bb bound states T(IS'), T(25') 
and T(35') with masses smaller than the threshold for BB production. The positions 
of these poles are known from the experimental masses of the T resonances, but 
the residues are unknown, containing the unphysical TBB or TBB* couplings 
The form factors V and Fi have in addition branch points at the threshold to of 
the BB production, below the beginning of the corresponding unitarity cut. If an 
estimate of the discontinuity across these cuts were available, the treatment of these 
subthreshold singularities in the present formalism could be done exactly (the 
method was applied recently in |r7| to the B D form factors). In what follows 



we shall resort to a pole approximation, keeping only the contribution of the T 
resonances situated below the thresholds tg and Q*, respectively. Using uib =5.279 
GeV, TUB* =5.324 GeV and the masses of the T resonances (Mx^ =9.460 GeV, 
= 10.023 GeV, = 10.355 and My^ =10.580 GeV) one can easily see that 
the form factor V{t) has only three poles below its unitarity threshold, much like 
F(t), while Fi(t) have four poles. Passing to the functions fi{z) according to (PD|) 



and using the conformal transformation ([T9|), we find that the functions fi{z) and 



/2(-2) have inside |z| < 1 three poles situated at the points 

2^ = -0.38, ;Z2 = -0.52,^3 = -0.67. (27) 

We neglected here the difference between and m^*, which is entirely justified 
as long as the singularities remain the same. As concerns the remaining functions 
fi,i > 3, they have four poles, with positions 

zl = -0.37 , z; = -0.49 , z; = -0.62 , ^ = -0.79, (28) 

obtained by using the conformal mapping (|T9|) with to replaced by and t by My-. 
The inequality (^^ has the form of an norm condition [10| involving several 



functions. We derive from it constraints on the functions fi and their derivatives at 
the origin z = 0, which corresponds through the conformal mapping to the point 
of zero recoil w = 1. If the functions were analytic, this would be very easily 
done, by applying standard techniques in the Hilbert space [|10|- However, as 
shown above, the functions have a finite number of poles, with known positions but 
unknown residua. The simplest treatment of this situation is based on the technique 
of Blaschke functions [l^ (the method was applied previously in |0]-@]). We define 
the following functions 

BW-At^^. B-{z)^Yll[^ (29) 

fc=i U ~ ^^k) U - zz,.) 

where we took into account that Zf^ and are real. 

As seen from ( p9D the functions B{z) and B*{z) have modulus equal to 1 on the 
boundary of the unit disk (i.e. for z = e*^). Therefore, we can insert the modulus 
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squared of the function B(9) (or B*{9)) in the integral appearing in (p^, without 
spoihng the inequahty or losing information. The relation (|2^ is thus equivalent to 

271^70 



^j2M(^)B^iO)f^m'de<l, (30) 



where we denoted 

B,iz) = B{z) (t = 1,2) ; Bi{z) = B*{z) {t = 3,6). (31) 

But the products Bi{z)fi{z) are functions analytic in \z\ < 1, the poles of the form 
factors fi being compensated by the zeros of the functions Bi{z). We can apply 
therefore the well-known results of interpolation theory for vector-valued analytic 



functions (see and references therein) to obtain from (|30D constraints on the 
form factors at points inside the analyticity domain. In particular, being interested 
in finding bounds on the charge radius (|^) which appear in (0), we shall apply an 
inequality of the Schur-Caratheodory type at z = 0: 

^ [<j).BJ,)\0) + {(j),BJ,y\0)] < 1 . (32) 

i=l 

It is important to note that up to now the form factors fi were treated as independent 
functions, without assuming that they coincide along the unitarity integrals. We use 
now heavy quark spin symmetry, which imply the relations (pT]). Then ( p2D can be 
written as an inequality for the charge radius 

5 5 

j2<pmBKo) + ^ [umm + 4>mBm - ^p^umm]' < i • (33) 

i=l 1=1 

The function /g does not contribute, due to the last condition in (^). The inequahty 
(p3D can be written as 

a{p^f - 26p2 + c < 0, (34) 

where 

a = 64^:5,^(0)0,2(0) 

i=l 

b = 8^^5,(0)0,(0) [0:(O)5,(O) + 0.(O)5:(O)] 

i=l 

5 5 

c = T.mo)Bm + <pmBmr + e Bfmm - 1- (35) 



i=l 1=1 



The quantities 0,(0), 0^(0), -8^(0) and -Bj'(O), entering the above coefficients, are 
calculable from the relations (^),(^) and ( 129|) and contain all the dynamical infor- 
mation in the problem. 
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3 Results and conclusions 



We discuss now the lower and upper bounds on the charge radius calculated 
from the above equation (|3^ . First we recall that the results previously reported 
in p[ and the second reference P] can be obtained by restricting the sums in the 



expressions (|35|) to a single term, i = 1. In the above works only the lowest order 
term 11'° in the expansion ( [T5| ) of IT' was retained and the contribution of the T poles 
in the relation (|T^ was dropped out. Also, for simplicity the choice mf, = tub for 
the mass of the h quark was made, and the value of which enters as a parameter 
in eq. ([131) was taken = 0. With these restrictions, eq. (0) gives the interval 
—4.5 < < 6.1 already reported in Keeping two terms {i = 1,2) in the sums 
appearing in (|35D , with the same numerical input as just described, we recover the 
interval —0.9 < < 2.60 obtained in |p. Finally, with all the five terms in the 
sums, i.e. by including all the form factors of the ground states B and B*, we 
obtain with the same input the range —0.35 < < 2.15. This result shows the 
improvement which can be obtained by fully exploiting spin symmetry in the ground 
state B doublet. 

As we mentioned, the above results were obtained with some simplifying assump- 
tions concerning the input. It is therefore of interest to perform the analysis with 
a more realistic input, according to the complete formulas given above. The main 
improvement is the QCD expression (p!5D of the polarization function up to three 
loops corrections. This expression depends on the scale p which appears in the MS 
coupling as{p) and in the coefficients C„ of the Taylor expansion (pTSl). We shall use 



in our analysis two scales, namely p = rrib, for which the coefficients Cn are [1^ 

Ci = 32.73, C2 = 33.24, = 29.61, Q = 26.94, (36) 

and p = 2mh, which gives 

Ci = 49.57, C2 = 43.31, C3 = 37.91, Q = 33.92. (37) 

We note that for the above choices of p the coefficients C„ do not depend on the 
specific value of rrib. Although the coefficients in ( |5BD and 0) are quite different, 
the final results, i.e. the bounds on p^, turn out to be practically the same. 

The expressions given in (15-18) were obtained using on shell renormalization. 



which means that is the pole mass. In the present work we shall treat this mass as 
a parameter in the reasonable range 4.7GeV — 5.GeV. For these values of rrib and the 
choices of p made above, the two-loops correction in the expansion (p!5| ) for = 
represents about 30% of the lowest order term, while the contribution of the three- 
loops diagrams is of about 10% (we used as(5.GeV) = 0.21 and ^^(lO.GeV) = 0.18 
). As we already pointed out, for heavy quarks one can extend the validity of 



the QCD perturbative expansion of the polarization function even at positive values 
of g^, below the threshold to. As an exemple, for g^ = 50 GeV^ the two-loops term 
represents a correction of about 45% of the one loop term, while the three loops 
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contribute in addition with approximately 20%. In the present formahsm better 
results, i.e. stronger bounds on p^, are obtained for larger q^. On the other hand, 
the increased contribution of the higher order QCD corrections for the polarization 
function prevents us taking too close to the hadronic singularities. We shall take 
in what follows q^ in the range — 50GeV^, noticing that the relative magnitude of 
the perturbative corrections does not dramatically change in this domain. 

We recall that much smaller values for the QCD perturbative corrections to the 
polarization function of heavy quarks were reported in [0 (see also ||13[)- The idea 
applied in these works was to express the pole mass in ([16|) and ( [17|) in terms 
of an euclidian mass defined to first order in as- This had the effect of reducing 
the procentual contribution of the two-loop correction, especially in the high order 
derivatives of the function n(q'^), of interest in the QCD sum rules for heavy quarks. 
The recent calculation of the polarization function up to three loops ||14| allowed 
us to use a more exact expression of 11', without resorting to the rather arbitary 



procedure adopted in |12 



The contribution of the T poles in the expression (JT^) was evaluated using the 
numerical values Vr^ = 1.34 keV, Tr^ = 0.56 keV and Tr, = 0.44 keV 111. The 
poles bring a positive contribution to the spectral function according to (|l3) and 
their inclusion improves the bounds in a significant way. 

In Fig.l we present the upper and lower bounds on the charge radius of the 
B meson elastic form factor, computed from (|3^), with the input described above, 
for mfy in the range 4.7GeV — S.GeV. As we mentioned, the two choices of the 
scale p adopted above give almost identical results. The solid curve corresponds 
to the choice = 0, the dashed one to = 50 GeV^. Taking larger values of 

we obtain much stronger bounds, but inconsistencies appear around 60.GeV^ 
(the pole contribution exceeds the QCD expression of n'(g^), signaling that a better 
estimation of the input is necessary). As seen from Fig.l, the predictions for the 
charge radius are rather sensitive to the value of the pole mass rrib, larger values of 
the mass leading to stronger bounds. 

The upper and lower bounds given in Fig.l represent the best results that can 
be derived, using a realistic input and fully exploiting heavy quark spin symmetry 
for the ground state B and B* mesons. We recall that the present derivation was 
possible by resorting to a a more powerful technique of analytic functions, which 
allowed the simultaneous treatment of several form factors as independent functions. 
The specific unitarity thresholds of the different form factors and their subthreshold 
singularities were correctly taken into account. Heavy quark spin symmetry was 
invoked finally by assuming that various form factors coincide near the zero recoil 
point, which is entirely legitimate. 

The technique applied in this paper can be easily generalized (see ||17[) to 



include higher derivatives of the form factors at the zero recoil point. In this way, for 
instance, quite strong correlations among the slope and the convexity of the elastic 
form factor F{t) can be derived. A second, more interesting generalization is to 
include in the unitarity sum the contribution of the excited states (B**) with orbital 
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momentum L = 1. By applying the techniques used in this work, it is possible to 
derive an inequality connecting the form factors of the ground states B and B* to 
the transition form factors between B** and the ground states. A new sum rule 
for these form factors, similar to the well-known inequalities of Bjorken and 
Voloshin PD[) will be reported in a future paper [ETI]. 
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Figure caption 

FIG.l: Upper and lower bounds on the charge radius of the elastic form factor of 
the B meson for various values of the pole mass mb. The solid line corresponds to 
= 0, the dashed one to = 50 GeV^. 



References 

E.de Rafael and J.Taron, Phys.Lett. B282, 215 (1992). 
[2] E.de Rafael and J.Taron, Phys.Rev.D50, 373 (1994). 
[3] I.Caprini, Z.Phys.C 61,651 (1994); Phys.Lett B339 187 (1994). 
[4] C.G.Boyd, B.Grinstein, R.F.Lebed, Phys.Lett. B353, 306 (1995). 
[5] N.Isgur and M.B.Wise, Phys.Lett.B232(1989) 113; B237(1990) 523. 
M. Neubert, Phys. Rep. 245, 259 (1994). 



[6] 
[7] 



[9 
[10 

[11 
[12 



N.N.Meiman, Sov.Phys.JETP 17, 830 (1971); V.Singh and A.K.Raina, 
Forts.der Physik 27,561 (1979); C.Bourrely, B.Machet, E.de Rafael, Nucl. Phys. 
B 189,157 (1981). 

P.Ball, H.G.Dosch, M.A.Shifman, Phys.Rev.D 47, 4077 (1993). 

I.Caprini Phys.Rev.D53, Nr.5 (March 1996). 

P.Duren, Theory of Spaces, New York: Academic Press, 1970. 

I.Caprini, I.Guiasu and E.E.Radescu, Phys.Rev.D 25 ,1808 (1982). 

M.A.Shifman, A.I. Vainshtein and V.I.Zakharov, Nucl.Phys.B147(1979) 385; 
B147 (1979) 448. 



11 



[13] L.J.Reinders, H.R. Rubinstein and S.Yazaki, Nucl.Phys. B186 109 (1981). 
[14] K.G.Chetyrkin, J.H.Kiihn, M.Steihauser, Phys.Lett. B371,93 (1996). 
[15] A.F.Falk, H. Georgi and B. Grinstein, Nucl.Phys. B343, 1 (1990). 
[16] I.Caprini, J.Phys.A:Math.Gen. 14, 1271 (1981). 

[17] Irinel Caprini and Matthias Neubert, preprint CERN-TH 95/225, accepted at 
Phys.Lett.B. 

[18] Review of Particle Properties, Phys.Rev.D50 1213 (1994). 

[19] J.D.Bjorken, in: Results and Perspectives in Particle Physics, Proceedings of 
the 4th Rencontres de Physique de la Vallee d'Aoste, La Thuille,Italy [SLAG 
Report, No. SLAG-PUB-5278, 1990, (unpubhshed)] 

[20] M.B.Voloshin, Phys.Rev.D46, 3062 (1992). 

[21] I.Gaprini, C.Macesanu and L.Micu, in preparation. 



12 



3 
2.5 

2 
1 .5 



0.5 



-0.5 

-1 



4.6 4.65 4.7 4.75 4.8 4.85 4.9 4.95 5 5.05 5.1 

m, (GeV) 



